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ABSTRACT 

The short pulse (SP) equation is a novel model equation describing the 
propagation of ultra-short optical pulses in nonlinear media. This article re- 
views some recent results about the SP equation. In particular, we focus 
our attention on its exact solutions. By using a newly developed method of 
solution, we derive multisoliton solutions as well as 1-and 2-phase periodic 
solutions and investigate their properties. 

1 INTRODUCTION 

In this article, we address the following short pulse (SP) model equation 

Uxt = U + ]:{U 3 ) XX , (1.1) 


where u = u(x,t) represents the magnitude of the electric field and subscripts 
x and t appended to u denote partial differentiation. The SP equation was 
proposed as a model nonlinear equation describing the propagation of ultra- 
short optical pulses in nonlinear media [1]. It is an alternative model equation 
to the cubic nonlinear Schrodinger (NLS) equation. The basic assumption 
in deriving the NLS equation is a slowly varying amplitude approximation. 
Hence, as discussed in the context of self-focusing of ultra-short pulses in non- 
linear media [2, 3], its validity would be violated if the pulse width becomes 
very short. A recent numerical analysis reveals that as the pulse length short- 
ens, the SP equation becomes a better approximation to the solution of the 
Maxwell equation when compared with the prediction of the nonlinear NLS 
equation [4]. Although the mathematical structure of the NLS equation has 
been studied extensively, only a few results are known for the SP equation. 
Here, we describe some recent results associated with the SP equation. In 
particular, we focus our attention on an exact method of solution, soliton and 
periodic solutions and their properties. 
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This article is organized as follows: In Sec. 2, we derive the SP equation 
starting with Maxwell equations for the electric and magnetic fields. In Sec. 
3, an exact method of solution is developed for the SP equation which trans- 
forms the SP equation to the integrable sine-Gordon (sG) equation through 
a hodograph transformation. In Sec. 4, the soliton solutions are constructed 
which include the multiloop soliton and multibreather solutions. Subsequently, 
the interaction process of solitons is described in detail. In Sec. 5, the exact 
method is applied to obtaining 1- and 2-phase periodic solutions. Some prop- 
erties of the solutions are discussed as well as their long-wave limit. In Sec. 6, 
an alternative method of solution is introduced which enables us to construct 
a more general class of periodic solutions. Then, the 1- and 2-phase solutions 
are exemplified. Section 7 is devoted to conclusion. 

2 SHORT PULSE EQUATION 
2.1 Basic equations 

The electric and magnetic fields E and H as well as the the electric and 
magnetic flux densities D and B are governed by the following set of equations 

<9B <9D 

div D = p, div B = 0, rot E = - — , rot H = j + — , (2.1) 

C6 C6 

where p and j are the electric charge and current densities, respectively. We 
consider the one-dimensional propagation of the wave so that we can put 

E = E 3 (x, t)e 3 , H = H 2 (x, t)e 2 , (2.2) 

where e 2 and e 3 are unit vectors perpendicular to the x axis. We also assume 
the following relations 

D = e E + P, B = p H, (2.3) 

where P is the induced electric@polarization, e is the vacuum permittivity 
and \i is the vacuum permeability. In view of (2.2), Eqs. (2.1) and the first 
equation of (2.3) are simplified to 

0H 2 _dD 3 0E 3 dH 2 

~dx~~^ ~dx~ ~ ^°~df' [2A) 

D 3 = e E 3 + P 3 , (2.5) 

respectively. Combining (2.4) and (2.5), we obtain the equation for E 3 . It 
reads in form 

E xx zE tt = P tt , (2.6) 

c z 
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where E = E 3 ,P = fioP^ and c 2 = (eo/Uo)^ 1 . The polarization P can be split 
into the linear part Pi in and the nonlinear part P n \ and it may be written in 
the form 

/oo 
X {1) (t-r)E(x,r)dT 
-OO 

/oo poo poo 
/ / X (3) (* - T lit - T 2 ,t - T 3 )E(x, T 1 )E(x,T 2 )E(x, T 3 )dTidT 2 dT 3 , 
-oo J —oo J —oo 

(2.7) 

where x^ an d X® are the susceptibilities. If we consider the propagation 
of light with the wavelength between 1600nm and 3000nm, then the Fourier 
transform x^ of x^ ls found to be well approximated by the relation x^ — 
X^ — X^ 1 ^ 2 [!]■ ^ follows from this and the relation u = 2ixc/\ that the linear 
equation for (2.6) written in Fourier transformed form becomes 

E xx + —^u 2 E - (2tt) 2 x 2 1) E = 0. (2.8) 

As for the nonlinear term in Eq. (2.6), we assume that the instantaneous 
contribution is dominant for the short and small amplitude pulses. Under this 
situation, we can set x^(t — r±,t — r 2 , t — r 3 ) = x^{t — Ti)S(t — r 2 )5(t — r 3 ) 
where X3 is a constant. If we introduce this relation into the nonlinear term 
on the right-hand side of (2.7), we obtain P nl = X3E 3 , which combined with 
(2.8), yields a single nonlinear wave equation for E 

E xx - \E tt = \e + X 3(E 3 )tt, (2.9) 

where c x = c/sjl + x^ and c 2 = l/(27ry // xl^). 

2.2 Perturbation analysis 

Equation (2.9) describes the interactions between the left and right moving 
pulses. Since the pulses are very short, the interaction between them would 
give rise to a higher-order effect on the evolution of the waves. Consequently, 
we may address only the right moving pulses, for instance. We use the multiple 
scale method to derive the approximate equation by expanding E as 

E(x,t) = eP o (0,X)+e 2 P 1 (0,X) + -- - , (2.10) 

where e is a small parameter which measures the shortness of the pulse relative 
to the time scale determined by the resonance, and and X are the scaled 
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variables defined by 

= — ^, X = ex. (2.11) 

If we introduce (2.10) with (2.11) into Eq. (2.9), we obtain, at the order 0(e), 
the following partial differential equation (PDE) for Eq: 

2 d 2 E n 1 „ d 2 E* 



ci d<pdX 6 



-£ + X3 —Ji. (2.12) 



After an appropriate change of the variables, we arrive at the normalized form 
of the SP equation (1.1). 

2.3 Remarks 

1. The SP equation has been derived for the first time in an attempt to con- 
struct integrable differential equations associated with pseudospherical sur- 
faces [5]. Schafer and Wayne rederived it starting from Maxwell's equations 
of electric field in the fiber as described in this section. See also a prior work 
due to Alterman and Rauch who discuss the breakdown of the slowly varying 
envelope approximation and perform an asymptotic analysis for a new type of 
nonlinear evolution equation [6]. 

2. The integrability of the SP equation has been established from various 
mathematical points of view [5, 7-10]. 

3. There exist several analogous equations to the SP equation which have been 
proven to be completely integrable. We write one of them in the form 

u xt = au+ ^(1 - (3)u 2 x - uu xx . (2.13) 

When (3 = 2, Eq. (2.13) becomes the short-wave model for the Camassa- 
Holm equation while when (3 = 4, it reduces to the the short-wave model for 
the Degasperis-Procesi equation and the Vakhnenko equation. The general 
multisoliton solutions of these equations have been obtained in parametric 
forms [11]. 

3 METHOD OF EXACT SOLUTION 

3.1 Reduction to the sine-Gordon equation 

Here we develop an analytical method for solving the SP equation which 
employs the hodograph transformation to reduce it to the completely inte- 
grable sG equation [12]. We first introduce the new dependent variable r 

r 2 = l+u 2 x , (3.1) 
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to transform the SP equation (1.1) into the form of conservation law 

1 



n = r j ■ (3.2) 
We then define the hodograph transformation (x, t) — > (y, r) by means of 

dy = rdx + -u 2 rdt, dr = dt, (3.3a) 

or equivalently 

^_ = r ^_ = _^_ + (3 36) 

dx dy dt dr 1 dy' 

In terms of the new variables y and r, (3.1) and (3.2) are recast into 

r 2 = 1 + r 2 u 2 y , (3.4) 

r T = r 2 uu y , (3.5) 
respectively. Furthermore, we define the variable <f> by 

u y = sm(f), = 0(y,r). (3.6) 

Inserting (3.6) into (3.4) gives 

- = cos0. (3.7) 
r 

It follows from (3.5)-(3.7) that 

u = <f) T . (3.8) 

Finally, if we substitute (3.8) into (3.6), we find that <fi obeys the following sG 
equation 

4> yT = sin0. (3.9) 

This form of the sG equation will be used to construct soliton solutions of the 
SP equation. For the periodic solutions, it is appropriate to introduce the two 
independent phase variables £ and r\ according to 

£ = ay + - + £o, (3.10a) 
a 
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T 

r] = ay \- r] , (3.106) 

a 

where a(^ 0), £o an d are arbitrary constants. In terms of the new variables, 
the sG equation (3.9) is transformed to 

0« - <f> m = sin 5 = </>(€> V)- ( 3 - n ) 

3.2 Parametric representation of the solution 

The solution u has a parametric representation given by (3.8). To be more 
specific 

u(y,T) = </> T . (3.12) 

To obtain the parametric representation of the coordinate x, we note from 
(3.3b) that the inverse mapping (y,r) — > (x, i) is governed by the system of 
linear PDE for x = x(y, r) 

x y = -, x T = -m 2 . (3.13) 

Since the integrability of the above system of equations is assured automati- 
cally by Eq. (3.5), we are able to integrate (3.13) immediately to obtain 

x(y,r) — J cos(f)dy + c, (3-14) 

where c is an integration constant. 

3.3 Criterion for the single-valued solutions 

As will be demonstrated later, most of the parametric solutions (3.12) 
and (3.14) become multivalued functions for both the soliton and periodic 
solutions. The single-valued functions are particularly useful in application to 
the real physical problem such as the propagation of nonlinear short pulses 
in an optical fiber. A criterion for single-valued functions may be obtained 
simply by requiring that u x exhibits no singularities. It follows from (3.3b), 
(3.6) and (3.7) that u x = tan0. Thus, if 

-7T<0<77, ( mod 7r )> (-V / 2+ 1 < tan^ < V2 - 1). (3.15) 

Zi £ J. 

then the parametric solutions (3.12) and (3.14) would become single-valued 
functions for all values of x and t. 

3.4 Remark 
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The reduction of the SP equation to the sG equation has also been estab- 
lished through the chain of transformations [8]. 



4 SOLITON SOLUTIONS 

4.1 Parametric representation of the iV-soliton solution 

In the context of the sG model, the soliton solutions are called kinks or 
breathers. These solutions are reduced from the soliton solutions by specifying 
the parameters such as the amplitude and the phase. Here we present the 
parametric representation for the iV-soliton solution of the SP equation [12]. 

The general iV-soliton solution of the sG equation can be written in a 
compact form as [13] 

* = 2iln£ (4.1) 

with 



/ 



n=o,i 



ex p 

l«=0,l 



N (fc + \ i) 



v^kijk 

l<j<k<N 



' N 

E 

" N 

.3=1 " l<j<k<N 



(4.2a) 



(4.26) 



Pj 



(j = l,2,...,N), 



(j,k = l,2,...,N;j^k). 



(4.2c) 



(4.2d) 



' Pj ~Pk \ 2 
K Pj+Pk) 

Here, p 3 - and £ j0 are arbitrary complex parameters satisfying the conditions 
Pj 7^ ±pk for j ^ k, i = \J— 1 and iV is an arbitrary positive integer. The 
notation X^=o i i m pli es the summation over all possible combination of \i\ = 
0, 1, fj,2 = 0, 1, hn = 0, 1. Note in (4.2c) that the variable r is replaced by t 
taking into account (3.3a). This convention will be used in the following. The 
r-functions / and /' play an essential role in constructing soliton solutions, as 
will be seen below. They satisfy the following system of bilinear equations 



ffyt fyft 



V 



A 



(4.3a) 



yt 



fyft 



J" 



i2 



n 



(4.36) 
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We obtain from (4.1) and (4.3) the important relation 

cos0 = l-2(lnfiV (4.4) 

Introducing (4.4) into (3.14) and integrating with respect to y yield the para- 
metric form of the coordinate x 

x(y,t)=y-2Qnf'f) t + c, (4.5) 

where c is an integration constant depending generally on t. It also follows 
from (3.12) and (4.1) that 

u(y,t) = 2i(\n^ , (4.6) 

which, combined with (4.6), gives the parametric representation of the N- 
soliton solution of the SP equation. To complete the solution, one must deter- 
mine the time dependence of c. To this end, we substitute (4.5) and (4.6) into 
the second equation of (3.13) and obtain the bilinear equation for / and /' 

fftt - V'tft + f'ftt = \c'{t)ff. (4.7) 

One can show that the r-functions / and /' from (4.2) vanish the left-hand 
side of (4.7). Consequently, d(t) = 0. Thus, in the case of soliton solutions, 
the constant c does not depend on t. In the following discussion, we consider 
real u so that we take /' = /* (complex conjugate of /). As will be seen 
later, this restriction imposes certain conditions on the parameters pj and 
£j (j — 1, 2, JV). The loop soliton and breather solutions of the SP equation 
are special classes of the general iV-soliton solutions given by (4.5) and (4.6). 

If we use (4.1) with /' = /* and the formula i ln(/*//) = 2 tan^Im //Re /), 
the criterion (3.15) for the single- valued solutions can be rewritten as 

-v / 2 + K^4<v / 2-l. (4.8) 
Re / 

4.2 Loop soliton solutions 
4.2.1 1-loop soliton solution 

Various solutions can be obtained for the SP equation by specifying the 
parameters pj and £jo(j = 1,2, ...,7V) in (4.2). The loop (antiloop) soliton 
solutions arise from the kink (antikink) solutions of the sG equation. Let m 
and iV-mbe the number of positive and negative pj, respectively. Then, 
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the corresponding soliton solution would describe the interaction of the m 
loop solitons and N — m antiloop solitons. Here we address the simplest 1- 
loop soliton solution which is of fundamental importance in discussing the 
properties of the general iV-loop soliton solution. In this case, (4.2) is written 

as 

/ = 1 + ie ?1 , 6 = PiV + —t + £io, (4.9) 

Pi 

2 . 5 

2 

1 . 5 

1 

. 5 



-6-4-2 2 4 6 
X 




Fig. 1: A typical profile of the 1-loop soliton solution. 

and /' = /*. Substituting (4.9) into (4.5) and (4.6), we obtain 

2 

u(y,t) — — sech£i, (4.10a) 
2 

x(y,t) = y - — tanh^ + d u (4.106) 

where d\ = c— (2/pi), pi > and ^i is a real constant. If we introduce a new 
variable X = x + c\t — x w with c\ = l/pf and x w = —^io/pi, then we can 
parameterize x(y,t) by a single variable ^. To be more specific, it reads 

X = — -— tanhfi + di. (4.11) 
Pi Pi 



It follows from (4.10a) and (4.11) that 

du sinh £i 



dX 2 - cosh 2 f ! ' 



(4.12) 
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We see from (4.12) that du/dX changes sign three times and goes infinity at 
£1 = ± cosh^ 1 2. Thus, the parametric solutoin exhibits singularities which has 
a form of single loop. The multi- valued feature of the solution is also confirmed 
by applying the criterion (4.8). In fact, it follows from (4.9) that Re / = 1 
and Im / = e ?1 and hence (4.8) cannot be satisfied for arbitrary values of 
y and t. Figure 1 shows a typical profile of the 1-loop soliton solution with 
the parameters p\ = l,di = 0. The loop soliton propagates to the left (i.e., 
negative x direction) at a constant velocity c±. If we define the amplitude A 1 
of the loop soliton by 2/p 1 (maximum value of u), then C\ = A\/A. Thus, the 
large loop soliton moves more rapidly than the small loop soliton, indicating 
the typical solitonic behavior. 

4.2.2 iV-loop soliton solution 

The general iV-loop soliton solution of the SP equation arises from (4.5) 
and (4.6) by taking the parameters pj(j = 1,2, ...,N) positive and £jo(j — 
1,2, ...,N) real. We first investigate the asymptotic behavior of the solution 
for large time and show that it is represented by a superposition of iV-loop 
solitons. The procedure for deriving the large time asymptptics can be per- 
formed straightforwardly by investigating the behavior of the r-functions given 
by (4.2). Hence, we omit the detail and describe only the result. To this end, 
we put Cj — 1/pj and order the magnitude of the velocity of each loop soliton 
as ci > c 2 > ... > c^. We observe the interaction of iV loop solitons in a 
moving frame with a constant velocity c n . We take the limit t — > — oo with 
the phase variable £ n being fixed. We then find the following asymptotic form 
of u and x: 

2 

u sech (f n + 5 { - ] ) , (4.13a) 

Pn 

2. . ^ . N 

x ~ 



y -l tanh (£ n + ,JH) - 4V--- + C. (4.136) 

Pn ^ +1 Pj Pn 

where 

N , s 2 

E ln(^^) • (4.13c) 

The corresponding asymptotic forms for t — > +oo are given by 

u ~ —sech (^ + 5 { n +) ) , (4.14a) 

Pn 

2 1 2 
x „ y - - tanh (£ n + 5^) - 4 V + c, (4.146) 

Pn ~(Pj Pn 
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with 

n-1 / \ 2 

4 +, = X>(^)- (414c) 

Let x c be the center position of the nth loop soliton in the (x, t) coordinate 
system. It then follows from (4.13) and (4.14) that 

s (-) N 1 

x c + c n t - x n0 ~ — 4 \-d n , (t -> -oo) (4.15a) 

P» Ai ^ 

n_1 1 

a; c + c n t - x n0 ~ — 4Y^ \- d n , (t -> +oo), (4.156) 

where x n0 = — £no/Pn and d n = c — 2/p n are phase constants. In view of the 
fact that all the loop solitons propagate to the left, we can define the phase 
shift of the nth loop soliton as 

A n = x c (t -> -oo) - x c (t -> +oo). (4.16) 

This quantity is evaluated using (4.13c), (4.14c) and (4.15) to give 

Pn [j^ \Pn+PjJ \Pn+Pj, 
(n-1 N \ 

+ 4 E - > (n = 1,2,..., AT). (4.17) 

Note that the first term on the right-hand side of (4.17) coincides with the 
formula for the phase shift arising from the interaction of N kinks of the sG 
equation. On the other hand, the second term arises due to the coordinate 
transformation (3.3). The latter changes the characteristics of the interaction 
process of loop solitons substantially when compared with those of the sG 
kinks. 

4.2.3 2-loop soliton solution 

The r-functions / and /' for the 2-loop soliton solution are written as 

f = 1 + ie ?1 + ie 6 - 7 e« 1+6 , (4. 18a) 

and /' = /* with 
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Fig. 2: The interaction of two loop solitons. 




Fig. 3: The profile of the 2-loop soliton solution. 

The parametric representation of the solution is then given by (4.5), (4.6) 
and (4.18). It reads 

t j.\ 2 V^ (Pi +Pz) cosh^iCosh^ 2 + (pi -P2) sinh^x sinh^ 2 , n N 
u{y,t) = — rn — ; r-r^- » (4.19a) 



P1P2 



cosh 2 ^1+7 sinh 2 ^ 



, 1 (pi -P2) sinh 2^i - t(pi + p 2 ) sinh2V;2 2(p 1 + p 2 ) t 

x \y^) — y + ,2 , , • 12 ; c ' 

P1P2 cosh tpi + 7 smh *0 2 P1P2 

(4.196) 
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where we have put 



^1 = ^(6-6), ?/>2 = ^i + 6) + ^ln7, (4.19c) 

for simplicity. The positive parameters pi and p 2 are assumed to satisfy the 
condition p 2 > pi- Figure 2 shows the interaction of two loop solitons with 
the parameters given by p\ = 0.5, p 2 = 1.0, c = £i = £20 — 0. 

Figure 3 shows the profile of the 2-loop soliton solution at t — —5 with 
the same parameters as those of Fig. 2. For the 2-loop soliton case, formulas 
(4.17) for the phase shift are written as 

A 1 = -llnf^^V-l, (4.20a) 

Pi \Pl+P2j P2 

A 2 = lln(^^y + ±. (4.206) 

P2 \Pl+P2j Pi 

Figure 4 plots P1A1 and piA 2 as a function of s(= P\jp2)- Thus, the large loop 
soliton always exhibits a positive phase shift whereas the small one exhibits a 
positive phase shift for < s < s c and a negative phase shift for s c < s < 1 
where s c is a solution of the transcendental equation A 2 = and is given by 
s c = 0.834. 




Fig. 4: The phase shift P1A1 and p\ A 2 as a function of s. The solid (broken) 
line represents the phase shift of the large (small) loop soliton. 
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4.2.4 Loop-antiloop soliton solution 

The solution representing the interaction of a loop soliton and an antiloop 
soliton arises if we choose p± > and P2 < with p\ < \p%\ in the 2-soliton 
r-function (4.2). The parametric solution takes exactly the same form as that 
of the 2-loop soliton solution (4.19). 




Fig. 5: The interaction of a loop soliton and an antiloop soliton. 



p l 




5 10 15 20 25 



x 

Fig. 6: The profile of the loop-antiloop soliton solution. 
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The formulas for the phase shift are given by 



Ax 


1 


In | 


(P\~P2 


Pi 




\P\ +P2 


A 2 


1 


In | 


(P\-P2 




Pi 




\P\ +P2 



+ — , (4.21a) 

P2 

+ — . (4.216) 
Pi 

Notice that the formula for Ai is altered when compared with (4.20a). Figure 
5 shows the interaction of a loop soliton and an antiloop soliton with the 
parameters given by p\ = 0.5, P2 = —1.0, c = £io = £20 = and Figure 6 shows 
the profile of the solution at t — — 5. 

4.3 Breather solutions 
4.3.1 1-breather solution 

The breather solution of the sG equation is the bound state of the kink and 
antikink solutions. Under certain condition, the breather solution is shown to 
yield a nonsingular oscillating pulse solution of the SP equation, which we 
shall term the breather solution as well. The 1-breather solution of the SP 
equation is derived if we put N = 2 in (4.2) and specify the parameters as 

Pi = a + ib, P2 = a — ib, (4.22a) 

60 = A + in, 60 = A - i/i, (4.226) 

where a and 6 are positive constants and A and \i are real constants. Then, 
(4.2) gives 

/ = 1 + ie ?1 + ie« + e« 1+ «, (4.23a) 

and /' = /* where £1 = 9 + i% with 

= a(y+ -^tj + A, (4.236) 
X = b(y-^—t) (4.23c) 



a 2 + 6 2 

Substituting (4.23) into (4.5) and (4.6), we obtain the following parametric 
representation of the solution 

4a6 6 sin vcosh (9 + In-) — a cos vsinh (8 + In-) 

u ^ * = -2—T2 h2 h 2 ( n a ' rrr 2 — r ' ( 424a ) 

a z + b z b 2 cosh [9 + In ^ J + a 2 cos 2 \ 
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x(y,t) = y 



lab asin2x + 6sinh2 (9 + In J) 4a 



a 2 + 6 2 6 2 cosh 2 [e + In J) + a 2 cos 2 x a 2 + 6 2 



+ a (4.246) 



Note that u has two different phase variables and The phase 9 charac- 
terizes the envelope of the breather whereas the phase \ governs the internal 
oscillation. In general, the solution (4.24) would exhibit singularities. Unlike 
the loop soliton solutions, the solution becomes a nonsingular function of x 
and t if we impose a condition for the parameters a and 6. To see this, we 
apply the criterion (4.8) to the r-function (4.2) to obtain 



V2 + 1 < 



a 



cosx 



b cosh (9 + In I) 



< V2-1. 



(4.25a) 



This inequality must be satisfied for any value of 9 and x- Since a > and 
b > 0, we see from (4.25) that the condition imposed on the parameters turns 
out to be 

0<a/6< v / 2-l. (4.256) 




Fig. 7: The profile of the nonsingular 1-breather solution. 

Figure 7 shows a profile of the 1-breather solution at t — with the pa- 
rameters a = 0.1,6 = 0.5, c = 80, A = /i = 0. In this example, a/6 = 0.2 
so that there appear no singularities as expected from the criterion (4.25b). 
For completeness, it will be instructive to present a singular breather solu- 
tion. Figure 8 shows an example of the singular solution with the parameters 
a = 0.4, 6 = 0.5, c = 20, A = /i = 0. Obviously, the criterion for the nonsingular 
solution is violated. 
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2 
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P 
-1 

-2 

~ 3 5 10 15 20 25 30 

x 

Fig. 8: The profile of the singular 1-breather solution. 

4.3.2 M-breather solution 

The general M-breather solution is constructed from the M-breather solu- 
tion of the sG equation (4.1) and (4.2) with N = 2M. Specifically, we set 

P2j-i = P* 2j = aj + ibj, aj > 0, bj > 0, (j = 1, 2, M), (4.26a) 

6,-1,0 = = \ + Vj, U = h 2, M), (4.266) 
and write the phase variables £2.7-1 and £ 2 j as 

= ^ + (J = 1, 2, M), (4.27a) 

£y = 0i-iXj> (J = 1,2,...,M), (4.276) 

with 

Oj = aj (y + Cj t) + Aj, (j = 1, 2, .., M), (4.27c) 

X 3 = bj(y - Cj t) + Hj, {3 = 1, 2, M), (4.27a 1 ) 

C ^ = ^T^' (J' = 1.2,...,M). (4.27e) 

The parametric solution (4.5) and (4.6) with (4.26) and (4.27) describes multi- 
ple collisions of M breathers provided that certain condition is imposed on the 
parameters aj and bj(j = 1,2,..., M). In the present M-breather case, the sim- 
ple inequality like (4.25b) is still difficult to obtain. However, as shown below, 
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the M-breather solution splits into M single breathers as t — > ±00. Hence, 
one can expect that the condition corresponding to (4.25b) would become 

M 

0< Y,1T <v / 2-l- (4.28) 
3=1 3 

It will be demonstrated later that the 2-breather solution exists whose param- 
eters indeed satisfy the inequality (4.28). 

Let us now investigate the structure of the M-breather solution by focusing 
on the asymptotic behavior for large time. To this end, we order the magnitude 
of the velocity of each breather as C\ > c<i > ... > cm- We take the limit 
t — > — 00 with 9 n being fixed. Then, we see that u and x have the leading- 
order asymptotics 

«<«.«) 



with 



An 



F n = b 2 n cosh 2 ^ n + + In ^ + a 2 cos 2 ( X „ + /3< _) ) , (4.29c) 

G„ = 6 n sin (xn + Pi^) cosh ^9 n 

-On cos (x„ + {3 { - ] ) sinh (0 n + a<"> + In — J , (4.29d) 

\ a n / 

# n = a„ sin 2 ( Xn + P { n ' ] ) + &n sinh 2 ^ n + + In ^ , (4.29e) 
where the real parameters ai, ^ and j3n are defined by the relation 



2M , s2 

E In(^V) =oir>+^ ) . (4.30a) 

The explicit expressions of ai ^ and /3 n ^ in terms of a,- and bj are calculated 
using (4.26a). They read 



, {(«» ~ S) 2 + (&n ~ kj) 2 }{K - S) 2 + (K + &j) 2 } . . 5()/ 

" {(a n + a j ) 2 + (6n + 6 J -) 2 H(«n + «i) 2 + (6n- W l ' J 
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/?<"> = 2 V f tan- 1 ^-^ 



_ tan- 1 - tan' 1 I . (4.30c) 

a n + a,- a n + dj I 

As t — > +00, -u and x take the same asymptotic forms as (4.29) with ai ^ 
and /5n ^ replaced by and respectively where 

2n-2 ✓ v 2 

The expressions of a4 + ^ and corresponding to (4.30b,c) follow if one re- 
places the sum YLjL n +i ^ Sj=i 111 (4-30b,c). Observing the asymptotic be- 
havior of the solution in the rest frame of reference, we see that it represents a 
superposition of M breathers, each has a form given by (4.24a). The effect of 
the interaction is the phase shift given by the sum of the quantities and 
(3^ which is caused by the pair wise collisions of M breathers and a term due 
to the coordinate transformation. 

The formula for the total phase shift will not be defined definitely since the 
solution takes the form of wave packets. If we consider the small-amplitude 
limit, however, the oscillating part and the envelope are shown to be separated 
completely so that one can obtain the formula like (4.17) for the iV-loop soliton 
solution. Actually, in the small-amplitude limit a n — > (n = 1,2, ...,M), 
expressions (4.29) and (4.30) are approximated by 

4 sin [Xn + /?n -) ) 

u(y, t) ~ ^ V J (4.32a) 

b n cosh (e n + a ( n ] + In y 

x(y, t)~y-^ tanh (e n + + In M - ^ + c, (4.326) 

V a nJ ° n 



" 8(o T 2 t + b])a 3 a n 

- fc2\2 ' {4.66(1) 



(b 2 - b 2 Y 

j=n+l 1 n 3 > 
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Let A n be the phase shift of the center position ( a point corresponding to the 
maximum amplitude) of the envelope for the nth breather. By performing the 
asymptotic analysis similar to that for the iV-loop soliton case, we find that 

A « = E L - E ■ (» = 1, 2, M). (4.34) 

J=7l+1 V n ^ i = l V " j> 

4.3.3 2-breather solution 

The solution describing the interaction of two breathers is parameterized 
by (4.26) and (4.27) with M — 2. Since the parametric solution has a lengthy 
expression, it is not appropriate to write it down here. Instead, we show 
the time evolution of the solution graphically and demonstrate its solitonic 
behavior. Fig. 9 depicts the profile of the two-breather solution for three 
different times (a)t = —40, (b)t = —5, (c)t = 35. The values of the parameters 
are chosen as a\ = 0.1, h\ = 0.5, a 2 = 0.16, 62 = 0.8, Ai = A2 = 0, Ui = U2 = 0. 
The velocity of the large breather is 3.85 while that of the small breather 
is 1.50 (see (4.27e)). Note in this example X^j=i( a i/^j) = 0-4 so that the 
inequality (4.28) is satisfied. For large negative time, the solution behaves like 
two independent breathers, each has a form given by (4.24) and propagates 
to the left. As time goes, both breathers merge and then they separate each 
other with leaving the original wave profiles. Apparently, the present example 
exhibits a typical feature common to the interaction of two-soliton solutions. 
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Fig. 9: The profile of the 2-breather solution for (a) t = —40, (b) t = —5, 
and (c) t = 35. 

4.4 Remark 

The 1- and 2-loop soliton solutions as well as the 1-breather solution have 
been obtained by different methods [14, 15]. 

5 PERIODIC SOLUTIONS 
5.1 1-phase solutions 
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The periodic solutions of the SP equation can be constructed by using an 
exact method of solution described in Sec. 3 [16]. Here, we deal with solutions 
which depend on the single variable 77. Then, the sG equation (3.11) reduces 
to an ODE for <fi 

0" = -sin0, (5.1) 

where the prime appended to (f> denotes the differentiation with respect to r\. 
There exist several particular solutions of Eq. (5.1). Among them, we look 
for solutions expressed by Jacobi's elliptic functions. Correspondingly, these 
solutions yield the 1-phase solutions of the SP equation. We investigate the 
properties of the solutions by the following two examples. 

5.1.1 Example 1 

The first example of the solution of Eq. (5.1) is given by Jacobi's sn 
function. Explicitly 

= ^sin^sn , (5.2) 

where the parameter k is the modulus of the elliptic function. Substituting 
(5.2) into (3.12) gives the parametric representation of u 

— I* (2.*), (5-3) 

where dn(w, k) is Jacobi's dn function. If we introduce the relation cos = 
1 — 2sn 2 (|, k) which is derived from (5.2) into (3.14) , we obtain 

x = y-2jsn 2 (J-,k^dy + c. (5.4) 

One can see that the integration constant depends on t whose time evolution 
is determined by the second equation of (3.13) with u given by (5.3). Indeed, 
using the identity 

jfe 2 S n 2 (|,fc)+dn 2 (|,jfc) =1, (5.5) 

we obtain an ODE for c, c'(t) = —2/(ka) 2 . This equation can be integrated 
immediately to give 

< t) = -W? t + lL (5 ' 6) 

where d is an integration constant. Substituting (5.5) and (5.6) into (5.4) and 
rearranging terms, we find a parametric representation of x 



I + W i - I » = 5{l ,2 - t, '" + 2£ (f t )} + ' 1 ' (5 ' 7 > 
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where x = —r] /a and E(u, k) is the elliptic integral of the second kind defined 
by [17] 



/" f T ll — kH 2 

dn 2 v dv = j y — — £2~dt, (t — sn v, r — sn u). (5.8) 

Thus, (5.3) and (5.7) give the parametric solution of the 1-phase solution. This 
solution becomes a multi-valued function. In fact, applying the criterion (3.15) 
to (5.2), we see that u x exhibits singularity when sn(rj/k,h) = ±l/y/2. To 
investigate the properties of the solution, it is convenient to take the amplitude 
and the modulus as independent parameters. The amplitude A of the wave 
may be defined by the relation A = (w max — u min )/2 where u max and u niin are 
the maximum and minimum values of u, respectively. In the present example, 
it follows from (5.3) that 

which enables us to express the parameter a in terms of A and k. The velocity 
V and the wavelength A of the wave are then given respectively by 

2-e (2 - e-)A' 

9 A 

A = 1 _ v / I — ^ {(2 - k 2 )K(k) - 2E(k)} , (5.11) 

where K(k) and E(k) are the complete elliptic integrals of the first and second 
kinds, respectively. In deriving (5.11), we have used (5.7) and the periodicity 
relation 

E(u + 2K(k), k) = E(u, k) + 2E(k). (5.12) 

The wavenumber K (which should not be confused with K(k)) and the angular 
frequency W are defined respectively by K = 2n/A and W = VK. The 
limiting forms of these wave parameters for k — > and k — > 1 are given 
respectively by 

K ~Tb> v ~ IF' w ~—' (k ^ 0) ' (5 ' 13a) 
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It follows from (5.13) that 



W ~ A 2 K (K -> 0), 
W ~ A 5 tf 3 /8 (K -> oo). 



(5.14a) 
(5.146) 




K 



Fig. 10: The dispersion relation W = W(K,A) with A = 1.0 for Example 1 
as a function of K. 




Fig. 11: A typical profile of the periodic solution for Example 1. 

The dispersion relation W = W(K,A) with A = 1.0 is plotted in Fig. 10 
as a function of K. Figure 11 illustrates the typical profile of the periodic 
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solution represented by (5.3) and (5.7) where the parameters are chosen as 
A = 1.0, k = 0.95, Xq = d = i] = and the time t is set to zero. In this 
example, one can see that the period A is 1.853. The figure represents a 
periodic loop traveling to the left at a constant velocity V = 2.320. 

When the wavelength of the periodic wave becomes very long, it degen- 
erates into a single loop soliton, as we shall now demonstrate. As seen from 
(5.13), the long- wave limit A — > oo (or K — > 0) is attained when k tends to 
1. Using the relations dn(w, 1) = sech u and E(u, 1) = tanh-u, the parametric 
solution represented by (5.3) and (5.7) reduces respectively to 

u = 2 A sech rj, (5.15a) 

x + AH + x = -Ar] + 2At&nh 7] + d, (5.156) 

with 7] = y/A — At + rjo. We see from this expression that the limiting solution 
is essentially the same as that of the 1-loop soliton solution given by (4.10a) 
and (4.11). See Fig. 1. 

5.1.2 Example 2 

The second example of the solution of Eq. (5.1) is given by Jacobi's dn 
function 

= — 2 cos -1 dn(r], k). (5.16) 
The parametric representation of the solution can be written in the form 

2k 

u = — cn(?7, k), (5.17a) 
a 



x - 4(1 - 2k 2 )t -x = -{-r] + 2E( V , k)} + d, (5.176) 



- 2k 2 )t - x = - 
a z a 

where cn(?7, k) is Jacobi's cn function. This solution is characterized by the 
wave parameters 

9k 

A=—, (5.18) 

(X 

K=^(1-2A; 2 ) = ^(1-2A; 2 ), (5.19) 



a 



2A 

A = — | - K(k) + 2E(k)\. (5.20) 

K 

Figure 12 shows the dispersion relation W = W(K, A) with A = 1.0. 

The dispersion curve has two branches depending on the value of k. The 
upper branch plotted by the solid line corresponds to the dispersion relation 
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for < k < k c whereas the lower one (broken line) represents the dispersion 
relation for k c < k < 1 where k c (= 0.9089) is a solution of the transcendental 
equation K{k) = 2E(k). Note that the wavelength A becomes zero when 
Jv — — k c . 
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2 
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-4 ; 

5 10 15 20 25 30 

K 

Fig. 12: The dispersion relation W = W(K, A) with A = 1.0 for Example 2 
as a function of K. 

The limiting forms of the wave parameters for both k — > and k — > 1 are 
given respectively by 

9k A 2 A 

K~^V~— W ~-(k^0). (5.21a) 

2-7T 1 A 2 IT A 

1 — k l — k 

We see from (5.21) that W ~ 1/fT (fT — > 0) for the upper branch and K ~ 
— A 2 K/A (K — * 0) for the lower branch. As K — > oo, both branches approach a 
straight line W = -V C K with V c = (2k 2 - l)A 2 /(Ak 2 ) ~ 0.1974A 2 . It is easy to 
see that the parametric solution (5.17) becomes a single- valued function when 
A; lies in the range < k < l/y/2. Note that the upper limit of this inequality 
coincides with the value of the modulus k for which the velocity given by (5.19) 
becomes zero. Figure 13 illustrates the profile of the nonsingular periodic 
solution at t — with the parameters A = 1.0, k = 0.65, Xq = d = i]q = 0. In 
this example, A = 3.027. It represents a periodic wavetrain traveling to the 
right at a constant velocity V = 0.0917. 
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Fig. 13: A typical profile of the nonsingular periodic solution for Example 2. 

If the parameter k lies in the range l/y/2 < k < 1, the solution exhibits 
singularities. Figure 14 illustrates the profile of the nonsingular periodic solu- 
tion at t — with the parameters A = 1.0, k = 0.8, x = d = i] = 0. In this 
example, A = 1.394. It represents a periodic wavetrain traveling to the left at 
a constant velocity V = 0.1094. 



3 




-2-10 1 2 



x 

Fig. 14: A typical profile of the singular periodic solution for Example 2. 

In conclusion, it will be worthwhile to consider the small amplitude limit of 
the solution. As suggested by the asymptotic relations (5.21a), the appropriate 
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limiting procedure is taken by the limit k — > while keeping the values of K, V 
and W finite. It turns out that the magnitude of the amplitude A is of order 
k. If we substitute the relations cn(i], 0) = cos rj and E(rj,0) = 77 into (5.17), 
the limiting form of the solution can be written as 

u = A cos ^ax - - -bj , (5.22) 

where b = a(x + d) is a phase constant. The dispersion relation of this linear 
wave is given by W = 1/K as is consistent with the asymptotics (5.21a). We 
also remark that (5.22) satisfies the linearized SP equation u xt = u. 

5.2 2-phase solutions 

5.2.1 Separation of variables 

The general iV-phase solution is now available for the sG equation. See 
[18], for instance. However, it will be difficult to perform the integral in (3.14) 
even for the 2-phase solution. An alternative approach for constructing the 
general iV-phase solutions will be discussed in Sec. 6. Here, we address the 
following specific form introduced by Lamb [19, 20] 



= 4 tan 



g(v) 



(5.23) 



We substitute (5.23) into the sG equation (3.11) and see that the variables £ 
and i] can be separated if / and g satisfy the following nonlinear ODEs 

f' 2 = - K f + nf + v, (5.24a) 

g ,2 = K g* + (n-l)g 2 -v, (5.246) 

where k, \i and v are arbitrary constants. For special choice of these param- 
eters, one can obtain solutions for / and g which are expressed in terms of 
elliptic functions. 

If we substitute (5.23) into (3.12), we immediately obtain the parametric 
representation of u 

«=^£ (5.25) 
a f 2 + g 2 

On the other hand, it follows from (5.23) by an elementary calculation using 
formulas of trigonometric functions that 

cos 0=1- * f * 9 * . (5.26) 
KJ +9) 
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The right-hand side of (5.26) can be modified in such a way that the integral 
in (3.14) can be performed analytically. To this end, we introduce the function 

" f 2 + g2 ' ^ 

where c\ and c 2 are constants to be determined later. Note that Y depends on 
the variables y and t through the relation (3.10). Now, we differentiate Y by 
y and use (5.24) to simplify the resultant expression. After some calculations, 
we obtain 



Y 



2K( Cl f + 3 Cl /V - 3c 2 /V - c 2 g G ) - 4C2/V 



+2(d + c 2 ) {-2/ ff /V + 2|i/V " -9 2 )}]- (5-28) 
We set ci + c 2 = and ci = — 2/a to reduce (5.28) in the form 

Y v = Mf 2 + 9 2 )--^^- (5-29) 

Comparing (5.26) and (5.29), we find that 

cos (f) = 1 + Y y - An{f + g 2 ). (5.30) 

Finally, we substitute (5.30) into (3.14) and take account of (5.27) with c\ = 
— c 2 = —2/a. Then, the integration with respect to y can be performed trivially 
to give the expression of x in terms of / and g 

x = y- -^4- -^1 U' 2 + 9 2 )dy + c. (5.31) 

The time dependence of c can be determined by the second equation of (3.13) 
with u and x given respectively by (5.25) and (5.31). It turns out that c'{t) = 
so that c = d(=const.). The expressions (5.25) and (5.31) provide the 
parametric representation of the the two-phase periodic solutions of the SP 
equation. 

We can obtain several 2-phase periodic solutions depending on the choice 
of the functions / and g. Here, we exemplify three solutions which reduce, in 
the long- wave limit, to breather solution (Example 1), 2- loop soliton solution 
(Example 2) and loop-antiloop soliton solution (Example 3). 

5.2.2 Example 1 
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The first example of / and g assumes the form 



f(0 = Acn(fc,k f ), g( v ) 



cn(Qi], kg) ' 



(5.32) 



where A, f3 and Q are positive parameters and kf and k g are moduli of the 
elliptic function. If we substitute (5.32) into (5.24), we can determine the 
parameters k, /i and v as well as kf, k g and Q in terms of A and (5. In 
particular, 

A ~ f l+ — 1 V (5.33a) 



1 + A 2 
A 2 



1 - 



1 + A 2 

n 2 = p 2 + 



P 2 (l + A 2 ))' 

1 \ 

Wil + A 2 ))' 

I- A 2 



1 + A 2 ' 



P 2 k 2 



(5.336) 
(5.33c) 
(5.33d) 



Note from (5.33) and the inequality < kf < 1 that the parameter f3 must be 
restricted by the condition 

- A < (3, (5.34) 

with arbitrary positive A. 

Now, the parametric representation of u follows from (3.12), (5.23) and 
(5.32). It reads 

_ AA -(3 sn(/3f , fc/)dn(/3£, k f )cn(Qr], k g ) + cn(/3£, k f )sn(Qr], k g )dn(Qr], k g ) 



(5.35) 



a A 2 cn 2 (/3£, fc / )cn 2 (^, fc p ) + 1 

The expression of x is dervied by substituting (5.32) into (5.31) as 



^2 cn(/3^, /c/)sn(f2?7, /c ff )dn(f2?7, /c 9 



a 



k f ) - fc} 2 # - {mv, k g ) - k> g 2 n V } 



+ d. (5.36) 
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Here we have used the following integral formulas for the Jacobi cn function 
in performing the integral in (5.51) 



/ 



J cn 2 (w, k)du = j-E(w, k) — A/ 2 w j , (5.37a) 

1 ( sn(u, k)dn(u, k) ,, 2 1 

-du = — 2 { ,\ - E(u, k) + k u \ , (5.376) 



cn 2 (u,k) k' 2 \ cn(u,k) 



where k! = y/T^W. Note that / (f 2 + g 2 )dy = a" 1 / + a' 1 / g 2 (rj)drj 

by virtue of (3.10). 

Let us now describe some properties of the parametric solution given by 
(5.35) and (5.36). In general, u is a multiply periodic function of x for fixed 
t. Under certain condition, however, it becomes a simply periodic function. 
To see this, we define the two parameters and L v by = 4K(kf)/j3 and 
L v = AK{kg)/Vt. In view of the periodicity of Jacobi's elliptic functions like 
sn(w + 4K(k), k) = sn(u, k), L^(L V ) is the period of u with respect to £(r)). In 
accordance with the value of A, there arise two possible cases for the period of 
u. When < A < A c (A c ~ 2.1797), one can show with use of (5.33) that the 
inequality L v < holds for arbitrary positive values of (3, and at fixed A, both 

and L T] are monotonically decreasing functions of (5 and vanish as (3 — > oo. 
If the ratio of both periods becomes a rational number, i.e., L^/L v = m n jm^ 
with (m^,m v ) = 1 and < m^, then u has a period L with respect to y 
given by 

L = —m^L^ = —m^L^. (5.38) 

(X (X 

With use of the relation (5.38), the period A with respect to x is determined 
from (5.36). Indeed, using the periodicity of the elliptic integral of the second 
kind E{u + 2mK(k),k) = E(u,k) + 2mE{k) (m : integer) as well as (5.33) 
and (5.38), the spatial period is found to be as 



A = L 



1-4/3 



1 



fE(k f ) kj E(k g ) 

\K(k f ) A*{l-kl)K{k a )^P{l + A*) 



(5.39) 



When A c < A, on the other hand, the equation = L v has a unique 
solution for j3 and the corresponding expression of the period is also given by 
(5.39). 

We remark that the solution presented here becomes a single-valued func- 
tion when the parameter lies in the rang < A < \[2 — 1. This follows from 
(4.8) and (5.32) with the aid of the inequality |cn(/3£, fcf)cn(fh/, k g )\ < 1. 
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Figure 15 depicts a profile of u at t — 0. The parameters chosen here are 
A = 0.2, = l,m v = 2, a = 1.0, d = £o — Vo — 0- Solving (5.38) for /3, one 
obtains (3 = 0.5832 so that Q = 1.124, fc/ = 0.3837, A; s = 0.0958, L = 11.21. 
Substituting these values into (5.39), the period A is found to be 10.37. 
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Fig. 15: A typical profile of the periodic solution for Example 1. 

We now consider the limiting profile of the periodic solution when the 
period A tends to infinity. To be more specific, we take the limits kf — > 1 
and kg — > 0. It then turns out from (5.33) that (3 — > A/y/1 + A 2 and f2 — > 
1/Vl + A 2 . The limiting value of (3 corresponds to the lower limit of the 
inequality (5.34). Using the relations 

sn(w, 0) = sin-u, cn(w, 0) = cosw, dn(w, 0) = 1, (5.40a) 

sn(w, 1) = tanhw, cn(w, 1) = sech u, dn(u, 1) = sech u, (5.406) 
E(u, 1) = tanhw, E{u, 0) = it, (5.40c) 
(5.35) and (5.36) reduces respectively to 



4AQ —A sinh 13 £ cos flri + cosh 3£ sin Qn 

u ~ s n o • (5.41a) 

cosh (3£ + A 2 cos 2 Qrj 



x^y ^ — L + a \ (5.416) 

a cosh (3£ + A 2 cos 2 Qry 
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This parametric solution is essentially the same as the breather solution al- 
ready given by (4.24). See Fig. 7. 

5.2.3 Example 2 

The second example is given by the following / and g 

m = A °J«§M ' M = _L ,5.42) 
cn(/3£, k f ) dn(ft?7, k g ) 



where 



*/ = 1 - A * + pif-A*) ' (5 - 43a) 

k ° = 1 - h + m=wy (5 - 436) 

n = (3A, (5.43c) 



K 



[3 2 (1 - kj) 



, n = f3\2 - kj), v = f3 2 A\ (5.43d) 



A 2 

The inequalities < kf < 1 and < k g < 1 impose the condition for f3 



I -A 2 ' 

with A in the range < A < 1. The exressions of u and x follows from (5.25), 
(5.31) and (5.42) with use of the formulas 

The resulting parametric solution reads in the form 

4A P dn(/3£, k f )dn(Qr], k g ) + k 2 £l sn(/3£, fc/)cn(/3£, k f )sn(Qrj, k g )cn(Qr), k g ] 

U = 

a A 2 sn 2 (^,A; / )dn 2 (fir/,A; 9 ) + cn 2 (^,fc / ) 

(5.46) 

_ 4(3 1 
^ ~ ^ ~ ~a A 2 sn 2 {(3£, k f )dn 2 {ttT], k g ) + cn 2 (/3£, kf) * 

x (A 2 dn 2 (^, kg) - l)sn(/3£, fc/)cn(/?£, %)dn(/3£, fc,) 
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+/^AW( / 3£, k f )sn(Qr], k g )cn(£lrj, k g )dn(£lr), k g ) 
+ M (-£(/?£, k f ) + AE(n V , k g )) + d. 

(X 



(5.47) 



Although the solution given above is a multiply periodic function, it has a 
single period if the condition 



L = lm € L € = j- a m v L v , 



(5.48) 



is satisfied. Using the relation (1 — kj)/k' g 2 = A 4 which follows from (5.43), 
the spatial period is found to be as 



A 



1 -A0- 



{ E(kf) 
\K(k f ) 



_ A 2 E(k 9 ) 



K(k g ) 



(5.49) 



Unlike Example 1, the solution always exhibit singularities as confirmed easily 
from (4.8) and (5.42). 

Figure 16 plots a profile of u at t — 5. The parameters are chosen as 
A = 0.2, = 2,m v — l,a — 1.0, d = £ — Vo = 0- I 11 this example, 
(3 = 1.027, Q = 0.2053, k f = 0.9998, k g = 0.8421, L = 20.35, A = 5.938. 




Fig. 16: A typical profile of the periodic solution for Example 2. 

In considering the limiting profiles, there arise two cases according to the 
inequality (5.44). The upper limit of the inequality for f3 is attained when 
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kg — ■> 0. In this limit, one has the limiting forms 



ft 



, fc f ~ 1 - A , / ~ A — 1 g~l. 



1-A 2,f ' ' ■'(•!!( K-kf) 

The expressions (5.46) and (5.47) then reduce respectively to 



(5.50) 



4 

u ~ - 



dn^M 



4 /M(l + A 2 ) dn^M 
aA 2 sn 2 (PZ,k f )+cn 2 (PZ,k f ) a dn 2 {^,k f ) + A 2 ' 

4(l + A 2 )sn(^,A; / )cn( / 9e,A; / )dn(^ ! A; / ) , 4/5 



(5.51) 



dn^A^ + A 2 



Using (3.10), (5.52) is modified in the form 



+^(-E((3^k f )+AQr 1 )+d, 

(5.52) 



I T/-J. _l n | ,02,2^ | 4 (1 + A 2 )sn(/^, fc/)cn(/%, fc/)dn(/%, M 

X+K ^o- a (l + 4^A J4 + a dn 2 (^,M + ^ 2 



A(3 



E(PZ,k f ) + d, 



with 



\ = \ +6A2 t A ! , *„ = -{(! + 4/3 2 A% - 4/3 2 ^ } 

a 




(5.53a) 



(5.536) 



Fig. 17: Periodic loop arising from the limit k g — > of the periodic solution 
depicted in Fig. 16. 
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Since the parametric solution (5.51) and (5.53) depends only on the single 
variable £, it becomes a 1-phase periodic function. Indeed, as shown in Fig. 
17, it represents a periodic train of loops propagating to the left at a constant 
velocity V. See also Fig. 1 which illustrates a typical profile of a 1-loop 
soliton solution. The maximum and minimum values of u are evaluated from 

(5.48) . They read as w max = 2(1 + A 2 )/[a{l - A 2 )] (at dn(/3£, kf) = A) and 
Umin = 4A/[a(l — A 2 )] (at dn(/3£, kf) = A 2 , 1). The spatial period A is given by 

(5.49) with L = K(kf)/(a/3). In this example, u max = 2.167, u min = 0.833, f3 = 
1.042, fi = 0.2083, k f = 0.9992, L = 4.442, A = 1.010, V = 1.347. 

The lower limit of (3 in (5.44) is realized when kf — > 1, fc 9 — > 1, which leads 
to the asymptotics 

fi ~ , f ~ A sinh/3£, (7 ~ cosh fir/. (5.54) 

V 1 — A 2 

Then, the parametric solution becomes 

4/3 A cosh /3£ cosh Q77 + A sinh /3£ sinh fir] 

it ~ o o 1 (5.55a) 

a A 2 sinh 2 /3£ + cosh 2 fir/ 



s ~ y - ^ ,r~, 2 " 2 - + d. (5.556) 



2/5 A 2 sinh 2f3£ - A sinh 2fir? 

a A 2 sinh 2 + cosh 2 fir/ 

This expression coincides with the parametric form of the 2-loop soliton solu- 
tion given by (4.19). See Fig. 3. 

5.2.4 Example 3 

The third example of / and g takes the form 

/ (fl = Xdn(«,t / ),^) = ^g, (5-56) 



where 



k * = 1 - i + ^(J-i) ' (5 - 57a) 

fi = |, (5.57c) 
K = ^, // = /3 2 (2 - fc 2 ), v = (3 2 A 2 (k 2 f - 1). (5.57d) 
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The inequalities < kf < 1 and < k g < 1 require that the parameter (3 
always must lie in the range 



A 2 -V 



(5.58) 



with A > 1. Using (5.25), (5.31) and (5.56), the parametric representation of 
the solution can be found to be as 



u = 



AA dn(/3£, kf)dn(Qi], k g ) + /3k 2 f sn(/3£, /c/)cn(/3£, kf)sn(Qr), k g )cn(flrj, k g 



x = y- 



A 2 dn 2 (/3£, k f )sn 2 {Qr], k g ) + cn 2 {Q V , k g 

4/3 1 

~a A 2 dn 2 ((3Z,k f )sn 2 (Qri, k g ) + cn 2 (^, k g ) 



(5.59a) 



x 



x — (1 - A 2 dn 2 ((3£, k f ))sn (Qrj, k g )cn (£lrj, k g )dn(ilr], k g 
-k 2 f A 2 sn^, k f )cn(PZ, fc,)dn(/3£, k f )sn 2 (Qr], k g ) 

(E(J3t, kf) - jE(n V , kSj + d. (5.596) 

As demonstrated readily by using (4.8) and (5.56), this solution always be- 
comes a multi- valued function. 

A spatial period of the above solution can be found if there exist in- 
tegers and satisfying the relation (5.48). Since in the present case 
K(kf)//3 < K(k g )/Q, one must impose the condition m v < m^, (m^,m v ) = 1. 
The expression of the spatial period is now given by 



A 



1-0 



fE(k f ) 1 E(k g ) 
\K(k f ) A 2 K(k g ) 



(5.60) 



Figure 18 plots the profile of u at t — 5. The parameters are chosen as 
A = 5,m^ = 2,m^ — l,a — 1.0, d = 12, £ = ?7o = 0. In this example, 
(3 = 1.027, Q = 0.2053, k f = 0.9998, k g = 0.8421, L = 20.35, A = 5.938. 

Last, we consider two limiting cases. When k g — > 0, (3 attains the upper 
limit of (5.58) and other parameters behave like 



A \/ A 4 — 1 

, k, ~ V A2 , / ~ A dn(/^, fc/), 9 ~ cot^. (5.61) 



A 2 - 1 
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The solution (5.59) reduces to 



u ~ 



AA ft dn(/%, k f ) + (3k 2 sn(/%, fc/)cn(/%, /c/jsinftr? cosftr/ 
" a A 2 dn 2 (/3£, A^sm 2 !^ + cos 2 ft?7 ' ' 



4/5 



a A 2 dn 2 (/3£, k f )sm 2 Qi] + cos 2 ft7? 



(1 - A 2 dn 2 (/3£, fey)) sin ft?? cos ftr? 



-k 2 f A 2 sn(^, k f )ca(M, k f )dn(J3t, k^sm 2 ^] - ^ (e(P£, k f ) - ifii^ + d. 

(5.626) 

In the case of = 2 and m v = 1, the relation (5.48) determines A uniquely as 
A = 6.553 so that (3 = 1.024, ft = 0.1562, = 0.9997, L = 20.11, A = 5.669. 
The solution (5.62) exhibits a profile similar to that depicted in Fig. 18. 




Fig. 18: A typical profile of the periodic solution for Example 3. 

The lower limit of j3 in (5.58) is established when kf — > 1 and k g — > 1. 
Consequently, one has 



ft 



=, / ~ A sech g ~ cosech flrj. 



Va^T' 

The solution then becomes 

4/3 cosh /3£ cosh ftr/ + A sinh /3£ sinh ftr? 



u ~ 



cosh 2 /3£ + A 2 sinh 2 ftr? 
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(5.63) 



(5.64a) 



2/? sinh2/%-Asinh2ta7 

rr ~ y s s \- «• 5.646 

a cosh 2 [3£ + A 2 sink 2 ttr] V ; 

It represents the interaction between a loop soliton and an antiloop soliton. 
See Fig. 3. 

5.3 Remarks 

1. An elementary method for obtaining 1-phase solutions is available which 
reduces the SP equation to a tractable ODE by assuming solution of traveling 
type [21]. 

2. The solutions (5.32), (5.42) and (5.56) have been derived in the context of 
a finite-length sG system [22]. See also [23, 24] for analogous works. 

6 ALTERNATIVE METHOD OF SOLUTION 
6.1 Bilinear transformation method 

The bilinear transformation method enables us to construct particular solu- 
tions of nonlinear evolution equations [25-28]. Although this method has been 
employed to obtain soliton solutions of the sG equation (see Sec. 4), it is appli- 
cable to periodic solutions as well. Indeed, Nakamura developed a systematic 
procedure for constructing periodic solutions of various types of soliton equa- 
tions [29, 30]. Here, we shall use his method to obtain periodic solutions of the 
sG equation. As already demonstrated in Sec. 4 for constructing soliton solu- 
tions, the r-functions play an essential role in the bilinear formalism. In the 
periodic problem, we introduce the same dependent variable transformation 
as (4.1) 

= 2iln^r. (6.1) 

Then, we can transform the sG equation (3.9) to the the following system of 
bilinear equations for the r-functions / and /' 

ffyt-fyft-\(f 2 -f' 2 ) = Xf 2 , (6.2a) 

rr yi -r y r t -\{f-f) = \r\ 

where A is a complex parameter to be determined later and the variable r has 
been replaced by the variable t by virtue of (3.3a). The parametric represen- 
tation of u follows immediately from (3.8) and (6.1) 

u(y,t) = 2i(\n f -^j . (6.3) 
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We then use (6.1) and (6.2) to derive the relation 

cos0 = l + 4A-2(ln/7V. (6.4) 

Introducing (6.4) into (3.14) and integrating with respect to y yield the para- 
metric representation of the coordinate x 

x(y, t) = (1 + AX)y - 2(ln f'f) t + c. (6.5) 

Comparing (6.5) with (4.6), one sees that a new parameter A comes in the 
periodic solution which would disappear in the long-wave (or soliton) limit. 
Thus, if we can solve the bilinear equations (6.2), then we can obtain solutions 
of the SP equation through the parametric representation (6.3) and (6.5). It 
should be remarked that unlike the soliton solutions, the constant c in (6.5) 
depends on t. This constant can be determined by using (4.7). 

6.2 Method of solution 

In accordance with Nakamura's procedure, we construct periodic solu- 
tions of the bilinear equations (6.2). To this end, we first introduce the N- 
dimensional theta function 

oo / N N \ 

0(z\t) = exp I 2ni^^njZj + ni n,jTj k rik j , (6.6) 

n\,n2,-.,n N =— oo \ j=l j,k=l / 

where z = (z±, Z2, z^) is an N- dimensional vector and r = (Tjk)i<j,k<N is 
an N x iV symmetric matrix. First, we seek solution of the bilinear equation 
(6.2a) in terms of the theta functions as 



(6.7a) 



where d = (1, 1, 1) is an iV-dimensional vector whose entries are all unity 
and Zj (j — 1,2, N) are phase variables defined by 



z 



■j 



kjy + Ujt + z j0 , (j = 1, 2, N). (6.7c) 

Here, kj,Uj and Zj are complex parameters. Substituting (6.7) into (6.2a), 
we find that the bilinear equation can be transformed to the form 

oo / n \ 

F(m,i, m,2, mjv)exp I 2ir\ S ^mjZj J =0, (6.8a) 

mi,m 2 ,...,mjv=-oo \ j=l / 
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where 



F(m 1 ,m 2 , ...,m N ) 



E 



rii,n2,...,npf=— oo 



N 



N 



-2tt 2 <^ ^2(2 nj - m 3 )k 3 \ \ J](2n, - m,) 



~[l + 4A- (-l)^=i^} 



x 



xexp 



A' 



N 



AT 



. (6.86) 



By shifting the sth summation index n s as n s + 1 in (6.8b), we see that 
F(m 1 ,m 2 , -,m N ) = -F(toi, ...,m s _i,m s - 2,m s+1 , ...,m N )x 



xexp 



JV 



Thus, if the relations 



\i=i 

F(m 1 ,m 2 , -,m N ) = 0, 



(6.9) 
(6.10) 



hold for all possible combinations of m 1 = 0, l,m 2 = 0, 1, ...,m N = 0, 1, then 
all F's become zero for arbitrary integer values of m 1 ,m 2 , ...,m N , implying 
that Eq. (6.8) holds identically. Consequently, the r-functions (6.7) satisfy 
the bilinear equation (6.2a). 

A similar analysis shows that the bilinear equation (6.2b) reduces to Eq. 
(6.8a) where the function F has the same form as (6.8b) except that the factor 
=1 rrij in the exponential function is replaced simply by — |i J2j=i m j- ^ 
turns out that the relations (6.10) assure that the r-functions (6.7) satisfy the 
bilinear equation (6.2b) as well. Thus, if we can determine parameters such 
that relations (6.10) are satisfied, then we obtain periodic solutions of the sG 
equation. We can regard (6.10) as a system of 2^ nonlinear equations for 
the unknown parameters u>j(j = 1, 2, N), Tjk(l < j < k < N) and A with 
given values of kj and Tjj(j = l,2,.,iV). The total number of unknowns is 
iV(iV - l)/2 + N + 1 = iV(iV + l)/2 + 1. For N = 1, 2, the total number of 
equations is equal to the total number of unknown parameters. Hence, we have 
1- and 2-phase solutions in terms of the theta functions. For iV > 3, on the 
other hand, the total number of equations always exceeds that of unknowns. 
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In this case, Eqs. (6.10) become an overdetermined system and we need a 
separate consideration as for the existence of the solution. 

6.3 1-phase solutions 

Here, we derive a 1-phase solution of the sG equation by means of the 
method described above. For N — 1, the relations (6.10) become 



F(m)= 



-2n 2 (2n - mfku - - {1 + 4A - (-1)™} 



x 



xcxp 



7T . 



7ri{(n — m) + n }r H — im =0, (m — 0, 1), 



(6.11) 



where we have put k = ki,u = u>i,t = Tu,m = m 1 ,n = n\ for simplicity. 
Explicitly, these read 

oo 

{^ 2 n 2 kcu + A)e 27rinV = 0, (6.12a) 

n=— oo 

oo 1 

{2vr 2 (2n - ifkuo + -(1 + 2A)}e" i{( "- 1)2+ri2}r = 0. (6.126) 

n=— oo 

We can rewrite (6.12) in terms of the following 1-dimensional theta functions 



e 1 (z\r) = -i (-l) n exp 



ni(2n + l)z + 7ri ( n + — ) r 



7ri(2n + l)z + vri ( n + - ) r 



02{z\t) = ^ eX P 
n=— oo 

9 3 (z\t) = exp [27rinz + 7rm 2 r] , 



n=— oo 
oo 



e 4 (z\r) = (-l) n exp \2-K\nz + vrinV] . 

n= — oo 

Before proceeding, it is convenient to introduce a new parameter q by 



(6.13a) 

(6.136) 
(6.13c) 
(6.13a 1 ) 



g = e" 



(6.14) 
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and write the above four theta functions as 



e j (z\T) = O j ( k Z,q), (j = 1,2,3,4). 



(6.15) 



Thus, the relation 9j(z\nr) = 9j(z,q n ) holds for any integer n. This notation 
will be used in the following. 

Now, using (6.15), Eqs. (6.12) can be recast into the following system of 
linear algebraic equation for the two unknowns kuj and A 



2#3(0,g 2 )Axc> -fl 3 (0,g 2 )A = 0, 



2^(0,g 2 )A^--# 2 (0,g 2 )(l + 2A) = 0, 



(6.16a) 
(6.166) 



where 9"(0,q 2 ) = d 2 9j(z,q 2 )/dz 2 \ z=0 , (j = 2,3). Solving this system, we 
obtain 

1 2 (0,q 2 )9 3 (0,q 2 ) 



ku = — 



4 9 2 (0,q 2 )9>>(0,q 2 ) - 9>>(0 : q 2 )9 3 (0,q 2 y 



A = -- 



9 2 (0,qyi(0,q 



2^(0,g 2 )^(0,g 2 )-^(0,g^ 3 (0,g 2 )' 



If we use the identity 



we can recast (6.17) to the compact expressions 

1 



ku = 



A = 



47r 2 fl 4 4 (0,g 2 )' 
9> 3 \0,q 2 ) 



2n 2 9 3 (0,q 2 )9i(0,q 2 y 
Now, the r-functions / and /' can be expressed in the form 



f = 9 ( Z + \ | T ) = £ 6XP 2nin ( Z + l) 



nin 2 T 



(6.17a) 
(6.176) 

(6.18) 



(6.19a) 
(6.196) 

(6.20a) 



T 



ex p 

n=— oo 



2 Tri// ( z — - ) -f Tiiir r 



(6.206) 
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with 



z = ky + ojt + zq. 



(6.20c) 



In order to obtain a real periodic solution, we introduce the new real quantities 
with tilde by 



i _ i ,r _ _ x 
z = -—z = - — {ky + Ut + Zq), 

Z7T Z7T 



(6.21a) 
(6.216) 



and put 

r = ib, (6 > 0), (6.21c) 
to assure the convergence of the series (6.20). Then, / and /' are rewritten as 



/= E 



exp 



n 



nn 2 b 



(6.22a) 



exp 



n 



7rn 6 



(6.226) 



with 



z = ky + out + z . (6.22c) 
In terms of the new parameters Co and k, the dispersion relation (6.19a) be- 

1 



comes 



(6.23) 



Obviously, /' = /* and A is real. Hence, the parametric solution given by 
(6.3) and (6.5) yields a real 1-phase periodic solution of the SP equation. For 
computing u from (6.3), we rewrite / in terms of the theta functions Q\ and 

64 as 



f = 04 



7T1 



i m 



7T1 



-,q 



(6.24a) 



Note that O^z/ni, q 4 ) and iQ^z/id, q A ) are real functions of The r-function 
/' is given by the complex conjugate of /. It reads 



7T1 



f' = e* 4,? 4 -i 4,? 



7T1 



(6.246) 
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It follows from (6.24a) and the definition of the sn function in terms of the 
theta functions that 

where 

T, = -i<? 3 W)*, «=H^|- (6-256) 
Furthermore, using the formula 

c \ • sn ( v ><> / n 2 ^1(°^ 4 ) /ro^ 

Sn If , K) = 1 : r, K = \J I — K 2 = -777 — , 6.26 

1 ^ cn(v,K'Y 01(0, q 4 Y 1 ' 

(6.25) becomes 

Imf 9 2 (0,q 4 )sn(9l(0,q 4 )z,K>) 



Re/ # 3 (0,g 4 )cn(fl 3 2 (0,g^X)' ^ 

The relation /' = /* makes it possible to write (6.3) as u = 4[tan _1 (Im//Re/)]t. 
Substitution of (6.27) into this expression yields u given by (5.51) with the 
identification among the parameters 

To derive the expression of x from (6.5), on the other hand, one needs the 
time dependence of c which can be determined from (6.24) and (4.7). After 
some calculations using the identities of the theta functions, we find 



c(t) = 



7T 2 



^4(0, q 4 



t + d, (6.29) 



where d is an arbitrary real constant. It can be demonstrated by substituting 
(6.24) and (6.29) into (6.5) that the expression (6.5) coincides with (5.52). In- 
deed, a straightforward calculation using (6.5), (6.20), (6.28) and some formu- 
las for the theta functions and Jacobi elliptic functions leads to the expression 
of x. We write it in the form 

4 (1 + A 2 ) sn(/%, k f )cn(ft, fc/)dn(^, k f ) 
V a dn 2 ((3^k f )+A 2 

+ t */> + ^wMk^) ") + 4Xy + c{t) - (6 - 30) 
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Note from (3.10a), (6.22c), (6.23) and (6.28) that ~z = (3£/9 3 (0,q 4 ) and from 
(3.10) that the last two terms on the right-hand side of (6.30) can be expressed 
in terms of £ and rj. With these facts in mind, we then use the formulas 



fl 3 '(0,g 2 ) _ 2fl 4 (0,gX(0,g 
91(0, qi) 



71 



£ 3 (0,g 2 )fl 4 (0,g 2 ) 2^ (0 ' g2) 



9 2 (0,q 4 )9 3 (0,q 4 



1 



* 2 W 



(6.31a) 



(6.316) 



and see that (6.30) coincides perfectly with the expression of x given by (5.52). 

Last, we consider the soliton limit. To this end, we first shift the phase 
constant z as z — > z + 7r6 and take the limit 6 — > oo (or o — > 0). In this 
limit, the theta functions 9 3 and #4 have the power series expansions 

9 3 (0,q) = l + 2g + 0(g 4 ), (6.32a) 

9 4 (0,q) = l-2g + 0(g 4 ). (6.326) 
Then, the asymptotic form of A from (6.17b) and that of uj from (6.23) become 

A ~ 0, uj ~ — . 

k 

The r-function from (6.22a) behaves like 



(6.33) 



/ ~ 1 + ie 5 , z 



r 1 
ky + -t 

k 



Zq. 



(6.34) 



This coincides with the r-function (4.9) for the 1-loop soliton solution. 

Another type of real 1-phase solutions can be constructed by a similar 
procedure to that described above. We list two of them for reference. If we 
replace z by z + \ and put r = \ + i6 in (6.20), the r-functions / and /' turn 
out to be 

00 r 

7T1 



/ = ^ exp 



n=— 00 

00 



27iinz 



-n 2 — irin 2 T 



7T1 



27rinz + — n" — mWr 



(6.35a) 



(6.356) 



where we have used the formula e nm ^ n+1 ^ — 1. In view of the formulas ^3(^|r + 
1) = 9^(z,r),9 4 (z\T + 1) = 9 3 (z,t), the relations (6.19) then become 



koj = 



4vr 2 # 3 4 (0,g 2 )' 



(6.36a) 
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A = 



27r^ 4 (0,g 2 )^(0,g 2 )- 



Using (6.35a), we find 



Im/ 
Re/ 



-A 



cn{w, K) 
dn(u>, k) ' 

Wg 4 )' 

By replacing z by iz in (6.35), we obatin 

Im/ A 

Re/ dn(w, k') ' 



where 

w = 2vr 2 ^ 2 (0,g 4 )(% + c;t + 2o), ^4 = 



« = A 2 , g = e 



-7T& 



(6.366) 



(6.37a) 



(6.376) 



(6.38) 



with k' = \/l — k 2 , where the expressions (6.36) remain the same forms. The 
expressions (6.37) and (6.38) give rise to the 1-phase solutions of the sG equa- 
tion. 

6.4 2-phase solutions 

In order to construct 2-phase solutions of the sG equation, we first solve the 
system of equations (6.10). Given values of k±, A; 2 , T\\ and T22, this system can 
be solved in principle for the unknown parameters uj\ , uji , T12 and A. Because 
of the transcendental nature of the system of equations, however, it is very 
difficult to obtain analytical solution. Under a few special situations in which 
the 2-dimensional theta function is expressed by a finite sum of 1-dimensional 
theta functions, the system can be solved algebraically. Indeed, there exist 
several examples to realize the situation mentioned above [31]. Among them, 
we consider the particularly important case 



Til = T 2 2- 



(6.39) 



It turns out that the 2-dimensional theta function (6.6) has the representation 
[31] 



om = - 



9, 



Z - + 2 



+ 9 3 (z + \t + )9 3 (z_\t_) 



where 



z± = -(z 1 ±z 2 ), 



(6.40a) 
(6.406) 
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t± = -(th ± r 12 ), (Imr ± > 0). 



(6.40c) 



It follows from (6.7) and (6.40) that 
1 



f = 0s 



i 

Z - + 2 



+ 9 3 



r = o 3 z+ + 



1 



+ 3 [z + - 



r + ^ 3 (^-|r_), (6.41a) 
r + ^ 3 (^-|r-). (6.416) 



Invoking the definition of the dn function, the solution <fi of the sG equation 
can be written in the form 



tan • 



lf-f 



4 if + f 

dn(w + + 5 + , k + ) y/kO — dn(w_, fc_) 
r k\ + dn(w + + 5 + , k + ) y / FT + dn(-u_, fc_) ' 



where 



\T±)Z±, 



k + = 



V 1 ~ k'±\ k' ± = 



|r ± ) : 



h = ^ 3 2 (0|r + ). 



(6.42a) 

(6.426) 
(6.42c) 

(6.42d) 



Thus, the above solution has a form similar to (5.23) in which the variables 
u + and U- are separated completely. 

The next step is to solve the system of equations (6.10) with N = 2. We 
observe under the condition (6.39) that 



nfru + 2n l n 2 r 12 + n\r 22 = (rii + n 2 ) 2 r + + (rii - n 2 ) 2 r_. 



(6.43) 



Using (6.43), Eqs. (6.10) can be solved analytically. The calculation involved 
is straightforward but somewhat tedious. Hence, we outline only the main 
steps. It follows from the equations F(l, 0) = and F(0, 1) = that 



kiUi = k 2 u 2 . 



(6.44) 



Substituting oo 2 from (6.44) into the equations F(0, 0) = and F{\,\) = 0, 
one obtains a homogeneous system of linear equations for uj\ and A. The 
solvability of this sytem yields the relation 

[e' 2 \o, q l)9 3 (o, q 8 + ) -e 2 (o, q lW(o, q l)][6 2 2 (o, q i) -el(o, q i)] 
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= a[0»(O, q*_)e 3 (0, q 8 _) - 9 2 (0, <£)0£(O, qt)} [£ 2 2 (0, £ ) - 61(0, £ )], (6.45a) 
where 

q± = e wiT ± = e^( T11±T12 )/ 2 , (6.456) 

Through a sequence of transformations using various formulas of the theta 
functions, one can simplify (6.45) to the form 

9l(0,ql) = a6 4 2 (0,q 2 _). (6.46) 

This is a transcendental equation which determines q + for given g_ and a. 
Once q + is obtained, the parameter r i2 is found from the relation 

which follows from (6.45b). The parameters and A are then determined from 
the equations F(0, 0) = and F(1,0) = as well as the relation (6.44), the 
explicit forms of which are not written down here. Thus, we have completed 
the construction of a 2-phase solution of the bilinear equation (6.2a). The 
corresponding solution of the SP equation can be obtained from (6.1) and 
(6.5). As in the case of 1-phase solutions discussed in Sec. 6.3, various 2- 
phase solutions of the SP equation arise by specifying the parameters ki,k 2 , 
and Tn(= T22). The explicit solutions are not presented here and will reported 
elsewhere. 

6.5 Remarks 

1. The starting point of our discussion is the condition (6.39) which makes it 
possible to perform all the calculations algebraically. The resulting solution of 
(j) has a separable form with respect to the independent variables u + and -u_ 
(see (6.42)). This expression should be compared with the separable solution 
introduced in (5.23). Thus, we can expect that the method developed here 
would reproduce all solutions constructed by a different method described in 
Sec. 5.2. 

2. Unless the condition (6.39) is imposed, we would be able to obtain a broader 
class of 2-phase solutions. The structure of solutions is worth studying in a 
future work. 

3. The general iV-phase solution of the sG equation has been constructed 
by means of the method of algebraic geometry [18]. It has a form given by 
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(6.1) in terms of the A-dimensional theta functions. However, whether the 
corresponding r-functions satisfy a system of bilinear equations (6.2) or not is 
an open problem to be resolved in a different context. On the other hand, our 
method first looks for the solution of (6.2) so that the expression of x follows 
immediately from (6.5). As already mentioned, however, the construction of 
the A-phase solutions with iV > 3 is a difficult task in the context of the 
bilinear formalism. Nevertheless, it is undoubtedly a challenging problem. 

7 CONCLUSION 

We have presented soliton and periodic solutions of the SP equation by 
means of a new method of solution. The most difficult technical problem 
in constructing solutions was how to integrate the PDE which governs the 
inverse mapping to the original coordinate system (see (3.13) and (3.14)). 
In the case of soliton solutions, the explicit form of the coordinate x was 
obtained in terms of the r-functions / and /' as shown by (4.5). In the 
case of periodic solutions, on the other hand, the special ansatz leads to the 
explicit form of x in terms of Jacobi's elliptic functions. Specifically, for 1-phase 
solutions, assuming the dependence of single variable, the sG equation becomes 
a tractable ODE (5.1) so that a number of special solutions are available. 
For 2-phase solutions, under the ansatz (5.23), we were able to reduce the 
representation of x to single integrals which are easily integrated (see (5.31)). 
It should be remarked, however, that the resulting solutions of the SP equation 
are special class of real 2-phase solutions. An alternative method employing 
the bilinear transformation method described in Sec. 6 enables us to construct 
a broader class of solutions than the solutions obtained in Sec. 5. At present, 
the most interesting issue will be the construction of the general A-phase 
solution which is reduced to the A-soliton solution (4.5) and (4.6) in the long- 
wave limit. 
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